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. Abstract 

Consider the single server queue with an infinite buffer and a FIFO discipline, either of type 
M/M/1 or Geom/Geom/1. Denote by A the arrival process and by s the services. Assume the 
stability condition to be satisfied. Denote by T) the departure process in equilibrium and by r the 
time spent by the customers at the very back of the queue. We prove that (2),r) has the same 
law as {A,s) which is an extension of the classical Burke Theorem. In fact, r can be viewed as 
the departures from a dual storage model. This duality between the two models also appears when 
studying the transient behavior of a tandem by means of the RSK algorithm: the first and last row 
of the resulting semi-standard Young tableau are respectively the last instant of departure in the 
queue and the total number of departures in the store. 

Keywords: Single server queue, storage model, Burke's theorem, non-colliding random walks, 
^ ■ tandem of queues, Robinson-Schensted-Knuth algorithm 

1 Introduction 

> 

I The main purpose of this paper is to clarify the interplay between two models of queueing theory. The 

first model is the very classical single server queue with an infinite buffer and a FIFO discipline. The 
second, less common but very natural, model can be described as a queue operating in slotted time with 

■ batch arrivals and services. It was studied for instance in |4j. The literature on discrete-time queues is 
' relatively less furnished than its continuous-time counterpart. The interest for discrete-time queues is 
. however gaining ground due to their relevance in ATM- like communication systems [5J ■ 

In this paper, to clearly distinguish between the two models, we choose to describe the second one with 
a different terminology and as a storage model. 

■ We prove first that the two models are linked in a very strong way. We set up an abstract model with 
^ I an ordered pair of input variables {A,s) and an ordered pair of output variables (T),r) = ^{A,s) = 

{^i{A, s),^2iA, s)). On the one hand, the queueing model corresponds to D = $1(^1,5) with A and 
s being respectively the arrivals and services and D the departures. On the other hand, the storage 
model corresponds to r = ^2iA, s), with s and A being respectively the supplies (arrivals) and requests 
(services) and r the departures. The interpretation of either r for the queue or T) for the store is much 
less natural. 

Then we assume that the random variables driving the dynamic are either exponentially or geometrically 
distributed, and we consider the models in equilibrium (under the stability condition). In this situation, 
it is well known that a Burke's type Theorem holds: the departures and the arrivals have the same law 
[71 inni m ■ This can be considered as one of the cornerstones of queueing theory, see for instance the 
books [51 [ini [31] for discussions and related materials. Here we prove a 'joint' version of the result: 
(D,r) has the same law as {A,s). This joint Burke Theorem is new, although it is similar in spirit to 
the results proved in [551 [H] for a variant: queues with unused services. Also, in the geometric case, 
we use an original method of proof based on the reversibility of a symmetrized version of the workload 
process (instead of the queue length process). As in [25l[T8], the joint Burke's Theorem can be used to 
obtain a representation theorem for the joint law of two independent random walks with exponential or 
geometric jumps conditioned on never colliding. 
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A second facet of the duality between queues and stores appears when considering the Robinson- 
Schensted-Knuth (RSK) algorithm. The RSK algorithm is an important object in representation theory 
and symmetric functions theory and it has been known for some time now that is has connections with 
queueing theory, although these connections have remained somewhat mysterious until now. In this pa- 
per we elucidate these connections completely via the above duality. The results we obtain, apart from 
being interesting in their own right, can be used as a starting point for the application of symmetric 
functions theory (and the associated analysis) to obtain interchangeability as well as asymptotic results 
for models which are constructed by putting the above queues and stores in series. These results are also 
applicable in the context of the corresponding interacting particle systems, where there is much interest 
in obtaining exactly solvable microscopic models for Burger- type equations. 

Let us now detail the results. Consider K queues (stores) in tandem: the departures from a queue (store) 
are the arrivals (supplies) at the next one. Initially, the network is empty except for an infinite number 
of customers (infinite supply) in the first queue (store). 

The particle systems associated with the above tandems are described in Section 15.21 One model is 
original and worth putting forward: the bus-stop model which is the zero-range model associated with 
stores in tandem. 

Here we assume that the variables driving the dynamic are N-valued without any further assumption. 
Building on ideas developed in [31 [53], we can study the transient evolution as follows. Consider the 
family of r.v.'s {u{i, j))i<i<js[^i<j<:K where u{i,j) is the service of the i-th customer at queue j, resp. 
the request at time slot i in store K + 1 — j. Apply the RSK algorithm, see [T71 [351 [M], to this family 
and let P be the resulting semi-standard Young tableau (here we do not consider the recording tableau 
Q). Let Ai > ••• > Xk > be the lengths of the successive rows of P. Classically [T7j, we have 
Ai = maxTren j)eir where 11 is the set of paths in going from (1, 1) to {N, K) and which are 

increasing and consist of adjacent points. Moreover, it is well known in queueing theory [TTl [5T1 [371 [35] 
that rnax^gn j)e7r instant of departure of customer N from queue K. Pasting the 

two results together gives the folklore observation that Ai = D. Here we complete the picture by proving 
that \k = R, where R is the total number of departures from the last store in the tandem up to time 
slot N. Again, this identity is proved in two steps by showing that \k — min^gj^ ^^-^ ^-^^^ m(j, j) = R, 

where H is a different set of paths in the lattice. To summarize, we obtain on the same Young tableau 
the total departures for the two tandem models. 

A short version without proofs of the paper appears in the proceedings of the conference Discrete Random 
Walks 03 [ig. 

2 Notations 

We work on a probability space (fl, P). The indicator function of an event i5 € 5" is denoted by 
1^;. We use the symbol to denote the equality in distribution of random variables. Depending on the 
context, \A\ is the cardinal of the set A or the length of the word A. We set N* = N\ {0}, = R+ \ {0}, 
and = a; V = max(a;, 0). We use the convention that 'Yl!i=j — Q when j > k. 

Below, a point process is a stochastic simple point process on R with an infinite number of positive and 
negative points. We identify a point process A with the random ordered sequence of its points: A — 
{An)n£Z with An < An+1 for all n. Observe that the numbering of the points is defined up to a translation 
in the indices. For any interval J, we define the counting random variable: ^1(1) = X^nsz ^{A„£i}- 
A marked point process is a couple {A,c) = {An,Cn)n(^7. where A = {An)ne'L is a point process and 
c = (c„)„gz is a sequence of r.v.'s valued in some state space. The mark Cn is associated to the point An 
of the point process. For precisions concerning point processes see [9]. 

Given a point process A — {An)nez, the reversed point process Ji{A) is the point process obtained 
by reversing the direction of time; i.e. 3l(yi) — {—A^n)nez- Given a marked point process {A,c) = 
{An,Cn)n(^i, the reversed marked point process is '5l{A,c) = (— A_„, c_„)„gz. 

Given a cadlag, i.e. right-continuous and left-limited, random process Y = (y(t))tgR valued in M, we 
define the reversed process "RoY = {"Ro Y(t))tesL as the cadlag modification of the process {Y {—t))ti^^. 
Denote by 3Nf+(y) and J^-iY) the point processes (with a possibly finite number of points) corresponding 
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respectively to the positive and negative jumps of Y, that is for any interval / of M, 

^+{Y){I) = l{Yiu)>Yiu-)}du, = ^ l[Yiu)<Yiu-)}du . (1) 

3 The Model 

Let A = {An)n&z be a point process and assume that Aq < < Ai. We define the ]R;!^-valued sequence 
of r.v's a = (a„)nGZ by a„ = An+i — An- Let s = (s„)„gz be another R!j_-valucd sequence of r.v's. The 
marked point process (^1, s) is the input of the model. 



Define the sequence of r.v.'s D — {Dn)n€Z by 



Dn = sup 

k<n 



Ak + (2) 



A priori the DnS are valued in R U {+oo}. Assume from now on that {A, s) is such that the Dn's are 
almost surely finite. They satisfy the recursive equations: 

Dn+i = max(i:)„, An+i) + s„+i . (3) 

Set T) = (-Dn)nez and set dn = Dn+i — Dn- Define an additional sequence of r.v.'s r = {rn)nei, valued 
in R;, by 

Tn = T[i\n{Dn,An+l) - An - (4) 

The marked point process (D, r) is the output of the model. By summing ^ and (H]), we get the following 
interesting relation between input and output variables: 

Tn + c?„ = a„ + s„+i . (5) 

In view of the future analysis, it is convenient to define the following auxiliary variables. Let w — (w„)„ez 
be the sequence of r.v.'s valued in R+ and defined by 

n-l 



Un=Dn-Sn-An^ SUp [V'(Si-ai) . (6) 
fc<n-l L ^ J 



i=k 

These r.v.'s satisfy the recursive equations: 

Wn+l = [Wn + Sn - On]^ • (7) 

Using the variables Wn, we can give alternative definitions of D„ and r„: 



< n, Dn = 



Wl 



n n 

Ai+Y^ s^] V max + ^ sj , (8) 



l<k<n 

i—l i—k 



Tn = min{u;„ + s„, a„} = s„ + w„ - w„+i . (9) 
At last, we define the cadlag random process Q ~ {Q{t))tes., valued in N, by 

Qit) = J2 Ma (10) 

nGZ 

Lemma 3.1. We have — A and ~ D. Furthermore, D^Qf^^^ < < -D_q(o)+i. 

Proof- Recall that Q{0) = X^nez 1{a„<o<Z5„} and that < < Ai. Hence, Q{0) = \k\ where k — 
sup{n eZ_\Dn<0}. The result follows. " □ 

We now interpret the variables defined above in two different contexts: a queueing model and a storage 
model. 
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3.1 The single-server queue 



A bi-infinite string of customers is served at a queueing facility with a single server. Each customer is 
characterized by an instant of arrival in the queue and a service demand. Customers are served upon 
arrival in the queue and in their order of arrival. Since there is a single server, a customer may have 
to wait in a buffer before the beginning of its service. Using Kendall's nomenclature, our model is a 
././1/oo/FIFO queue. 

W(t) . 




Figure 1: The dual variables (s„)„ and (r„)„ 

The customers are numbered by Z according to their order of arrival in the queue (customer 1 being the 
first one to arrive strictly after instant 0). Let An be the instant of Arrival of customer n and s„ its 
Service time. Then the variables defined in (|^- (fTU)) have the following interpretations: 

• Dn is the instant of departure of customer n from the queue, after completion of its service; 

• Wn is the waiting time of customer n in the buffer between its arrival and the beginning of its 
service; 

• Q(t) is the number of customers in the queue at instant t (either in the buffer or in service); 
Q = {Q{t))t is called the queue-length process; 

• r„ is the time spent by customer n at the very back of the queue. 

The variables (r„)„ are less classical in queueing theory although they have already been considered |26| . 
They should be viewed as being dual to the services (sn)n as illustrated in Figure [TJ On the upper part 
of the figure, we have represented the workload process {W{t))t, where W{t) is the waiting time of a 
virtual customer arriving at instant t (see (|12|) for the formal definition). 

3.2 The storage model 

Some product P is supplied, sold and stocked in a store in the following way. Events occur at integer- 
valued epochs, called slots. At each slot, an amount of P is supplied and an amount of P is asked for by 
potential buyers. The rule is to meet all the demand, if possible. The demand of a given slot which is not 
met is lost. The supply of a given slot which is not sold is not lost and is stocked for future consideration. 

Let Sn be the amount of P Supplied at slot n + 1, and let a„ be the amount of P Asked for at the same 
slot. The variables in ((2|l- (fT0|) can be interpreted in this context: 

• Wn is the level of the stock at the end of slot n. It evolves according to ([7|); 
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• r„ is the demand met at slot n + see equation ([9]); it is the amomit of P departing at slot n+1; 
The variables {Dn)n a-nd the function Q do not have a natural interpretation in this model. 

The evolution of the store is summarized in Figure^ The indices may seem weird (s„,a„,r„, for time 
slot n + 1). They were chosen that way to get better looking formulas in fj^l 

Time slot n Time slot n+1 

Sn_i a||_| rj|_| Sj^ 

Figure 2: The storage model 

It is important to remark that while the equations driving the single server queue and the storage model 
are exactly the same, it is not the same variables that make sense in the two models. The important 
variables are the ones corresponding to the departures from the system. The departures are coded in 
the variables (Dn)™ for the single server queue and in the variables (r„)„ for the storage model. On 
the other hand, interpreting the variables (r„)„ in the single server queue or the variables {Dn)n in the 
storage model is not so immediate. 

To summarize, the output variables and (r„)„ are equally relevant, but in different modelling 

contexts. 

Observe that it would be possible to describe the above storage model with a "queueing" terminology (a 
queue with slotted time, batch arrivals (s„) and batch services (a„)). It is the description used in [4j for 
instance. We have avoided it on purpose to clearly separate the models of §3.11 and §3.21 and therefore 
minimize the possibility of confusion. 

4 Equilibrium Behavior: Around Burke's Output Theorem 

We consider the exponential or geometric version of our model in equilibrium. We prove a Burke's type 
result: {T>,r) ^ {A, s). The relevant result for the sequence of departures in one of the two models will 
then follow by forgetting one of the two variables. In the exponential case, our proof uses the queue- 
length process Q following [30] . Our method of proof in the geometric case is original and based on the 
zigzag process, a symmetrized version of the workload process. A discussion of the literature and of the 
increment of the present version is carried out in §4.31 

4.1 Output theorem in the exponential case 

Let A be an homogeneous Poisson process of intensity A £ R*^. We set A — (A„)„gz with Aq < < Ai. 
Recall that a„ = A„+i — An- Then (a„)„>i is a sequence of i.i.d. r.v.'s with exponential distribution 
of parameter A. Let s — {sn)n<£Z be a sequence of i.i.d. random variables, independent of A, with 
exponential distribution of parameter /i. We assume that A < /i. 

Consider the marked point process (^L, s) as being the input of the model of Section [3l The sequence 
{wn)n is a random walk valued in R_|_ with a barrier at 0. Under the stability condition A < /z, this 
random walk has a negative drift. It implies that the random variables D„ defined in ([2]) are indeed 
almost surely finite. 

We now prove Theorem 14.11 Figure [2] provides a visual illustration of the result. The principle of the 
proof is the same as in Reich [30j . However, the adaptation is not completely immediate. Some care is 
necessary in dealing with the indices, to make sure that the variable r„ is the mark of the point Z?„. 

Theorem 4.1. The marked point process (33, r) has the same law as the marked point process {A, s). 
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Proof. Knowing Q, one can recover the input of the model: (An, Sn)n = '^{Q)- We are going to prove 
that 

ipo'3l{Q) = (-£'-„+Q(o)+i,r_„+Q(o)+i)„ =3?('D,r) 

Let T„ be the operator mapping a point process to its n-th point, with the convention that To(-) < < 
Ti(-). For instance, An = 'Jn(yi) = T„ o ?\f+((5). Taking into account the shift in the indices and using 
Lemma I3TT1 we get that T„(2)) = T„ o yi_{Q) — Dn-Q(o)- So we have 

= %i+Q{0) ° ^-(Q) - niax(T„ o K+(Q), T„+q(o)_i o 'N-iQ)) . 
Set Sn ~ §n{Q)- Now let us apply the operators T„ o 3\f+ and §„ to the reversed process 

T„ O O Ji{Q) = -D_n-Q(0) + 1 

Furthermore Tn+JiiQW)-! oT^-o ^{Q) = -^-«-q(o)+2- Thus, 

S„o3l((5) = -A_„_Q(o)+i - max( -Z3_„_Q(o)+i, -A_n-Q{o)+2) 

= -A_n+Q(Q) + 1 + niin( D_„+Q(o) + i, ^_„+g(o)+2 ) 

= r-n+Q{0) + l ■ (11) 

Hence we obtain that (p o = 31(1), r). 

The process Q is a stationary birth-and-death process, hence reversible: ^{Q) ^ Q- It implies that 
{A,s) = f{Q) ^ (f o 'Ji{Q) = 3^(D,r). Hence, Jl{D,r) is an homogeneous Poisson process marked with 
an i.i.d. sequence of r.v.'s. So, 3?(D,r) ~ (I',r), which concludes the proof. □ 

Corollary 4.2. In the queueing model, the departure process D is a Poisson process of intensity A. In 
the storage model, the sequence (r„)„ of the amounts of product P departing at successive slots, is a 
sequence of i.i.d. exponential r.v.'s of parameter ^. 

4.2 Output theorem in the geometric case 

Let be a Bernoulli point process of parameter p S (0, 1), that is: all the points are integer valued, 
there is a point at a given integer with probability p, and the presence of points at different integers are 
independent. As before, set A = with Aq < < Ai and a„ = An+i — An. Then the sequence 

(an)n>i is a sequence of i.i.d. geometric r.v.'s with parameter p (Vfc G N*, P{ai = fc} = (1 — p)''~^p). 
Let {sn)n be a sequence of i.i.d. geometric r.v.'s with parameter q e (0, 1), independent of .A. We assume 
that p < q (stability condition). 

Let the marked point process {A,s) be the input of the model of Section [3l As in HA.ll the model is 
stable and the output (©jr) is a marked point process. 

We prove a Burke-type Theorem using an original approach based on a symmetrization of the workload 
process. 

Define the (random) cadlag maps /„, g„ : f2 x K ^ M, by 

'O ift<An fO ift<An 

fn{t)^{Dn-t iiAn<t<Dn, gn{t)=\t-An liAn<t<D,, 



Set 



iit> Dn iit> D. 



W{t) = V /„(t), W{t) = V gn{t) . (12) 

Since Dn — An + w„ + s„, we have W{An ) = w„. For the queueing model, W{t) is the total amount of 
service remaining to be done by the server at instant t, and W = {W{t))t is called the workload process. 
The process W = {W{t))t is the dual of the process W , see Figure [T] or Figure [3] for an illustration. 
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Looking at W and W together is a way to symmetrize the workload process. Consider the ordered pair 
of processes P = {W, W). Define 

51{P) ^ {5l{W),5i{W)) 

(observe that the first element of the pair is now Jl{W), as opposed to W for the pair P). Obviously, we 
can recover the input of the model knowing P: {A, s) = 'I'(-P). Also it is clear (see Figure [T] or Figure [3]) 
that: 3?(D,r) = 4'o3?(P). 

The crux of the argument in proving a Burke-type Theorem is now the next Lemma. 

Lemma 4.3. The ordered pair of the workload process and its dual is reversible, that is: [W, W) ~ 
{Ji(W),5liW)). 

Proof. Clearly, W and W are (time) stationary and ergodic processes. An idle period is a maximal 
non-empty time interval during which W{t) = (or equivalently, W(t) = 0). The busy periods are the 
time intervals between the idle periods. Hence R is partitioned by an alternating sequence of busy and 
idle periods. 

Let (in)nez and (6„)„gz be the lengths of the successive idle periods and busy periods respectively. 
Clearly, the sequences («„)« and (6„)n are independent and i.i.d. Moreover, it follows from the memoryless 
property of the geometric distribution that i„ is a geometric distribution of parameter p. So, to prove 
the reversibility of {W, W), it is enough to prove the reversibility within a single busy period. 

To do so, the easiest is to introduce an auxiliary process. Consider a busy period which is assumed to 
consist, for simplicity, of the k customers numbered from 1 up to k. Set for n ^ {1, . . . k}: 

Cn = Bn + s„, Bn+i = C„ + a„ . (13) 

This defines (i3„)n and (C„)„ up to a translation. Define 

^ ' ^' \[Z{Cn) - {t - Cn)]^ forte [C7„,S„+i) ^ ^ 

On an interval of type [Bn,Cn), we have dZ/dt = 1 and on an interval of type [Cn, Bn+i), we have 
dZ/dt={-l)lz>o- 

Now, on a busy period, the process Z is in bijcction with {W, W), and the time-reversed process 3{{Z) is 
in bijection with {31{W) , Jl{W)) . This is illustrated in Figure[3l where we have represented a trajectory 
of {W, W) and the corresponding trajectory of Z. The process Z is obtained by gluing together the slices 
of W and W alternatively. 




Figure 3: The workload process, the dual workload process, and the zigzag process 

Hence, proving the reversibility of (W, W) on a busy period, is equivalent to proving the reversibility of 
Z. We may call Z the zigzag process. This process was considered (with a different motivation) in [12j 
and in particular, it was proved there that the zigzag process is reversible (although not Markovian). 
We recall the argument for completeness. 
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Identify Z with the finite sequence of the lengths of its intervals of increase and decrease. Let 
(/i,...,/2fe) G (N*)^*^, with X^i = Si'2i+i = L, be a possible value for Z. (The total length of 
the periods of increase has to be equal to the total length of the periods of decrease). Recall that the 
variables (sn)n and (a„)„ are geometrically distributed of respective parameters q and p. We get: 

p{z - (Zi, . . . , hk)) = (1 - ■ ■ ■ (1 - 9)'^'="' 

- (l-p)^-V(l-g)^-V-^ 

Hence, given L and fc, all the different possible trajectories for Z arc equiprobable. It implies that the 
same is true for the time- reversed process This completes the proof. □ 

We are now ready to prove the following result: 

Theorem 4.4. The marked •point process (Ti,r) has the same law as the marked point process {A,s). 

Proof. The proof is similar to the one of Theorem 14 . 1 1 wit h {W, W) playing the role of Q. 
Recall that {A, s) = ^{W,W) and 3?(D,r) = o Ji{W,W). According to LemmaEJl we have 

{A,s) = ^'(W^,TF) ~ *o3?(VF,TF) = Jl{V,r) . 

We conclude that (D, r) ^ Jl{A, s) ~ {A, s), where the last equivalence comes from the fact that {A, s) 
is a Bernouilli process with i.i.d. marks, hence is reversible. □ 

4.3 Comments on the different proofs of Burke Theorem 

Reflecting on the above, there are three different ways to prove Burke Theorem, be it in the exponential 
or geometric case. The first way is by using analytic methods, the second is by using the reversibility of 
the queue length process Q, and the third is by using the reversibility of the zigzag process Z. 

The original proof of Burke is for the exponential model using analytic methods [7] . For the geometric 
model, an analytic proof was given by Azizoglu and Bedekar [3]. For the exponential model, the idea of 
using the reversibility of Q to get the result is due to Reich [3D] . This proof has become a cornerstone of 
queueing theory, it has been extended to various contexts and has given birth to the concept of product 
form networks O [15]. Reich's proof does not translate directly to the geometric model. Of course, 
{Q{n))nei, is a reversible birth-and-death Markov chain. However, a difficulty arises: It is not possible to 
reconstruct A and s from Q. Indeed, on the event {Q{n — 1) = Q{n) > 0}, two cases may occur: there is 
either no departure and no arrival at instant n, or one departure and one arrival; and it is not possible 
to distinguish between them knowing only Q. One feasible solution is to add an auxiliary sequence that 
contains the lacking information but the details become quite intricate. This program has been carried 
out in [13], see also [TH] Theorem 4.1], for a variant: the geometric model with unused services. The 
above idea of using the zigzag process to prove Burke Theorem is original. Observe that the zigzag 
process is also clearly reversible in the exponential model. Hence the proof used for Theorem 14.41 could 
also be used to get Theorem 14. II 

In the original references [T] [3^ and in all the classical textbooks presenting the result [5] [HI [31] , the 
version proved is: A ^ T>: ("Poisson Input Poisson Output"). In [3], or [T3| for the variant with unused 
services, the result proved is s ~ r. The complete version (A, s) ~ (D, r) appears first in [25, Theorem 
3] for the exponential model with unused services. This is extended to the geometric model with unused 
services in [18]. Brownian analogues are proved in [131 124j . 

4.4 Non-coUiding random walks 

Following the lines of thought in [251 [HI [E] , it is possible to use Theorems 14.11 and 14.41 to get represen- 
tation results for non-colliding random walks. 

We start by stating a preliminary lemma. 

Lemma 4.5. The Lindley's Equation ^ admits a backward analogue: 

Wn = {wn+i + r„ - d„_i)+ . (15) 
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Proof. Recall Equation ([9]): r„ = s„ + w„ — Wn+i- Using ([5]), we get 

rfn = On + S„+l - r„ = a„ - W„ + Wn+1 - Sn + S„+i . (16) 

Therefore, 

Wn+l + rn - dn-1 = W„+i + S„ + W„ - W„+i - a„_i ~ w„ + w„_i - S„ + S„_i 
= Wn-l + Sn-1 - a„_i . 

Taking the maximum with on both sides of the equaHty and using Lindley's equation ((7|), we obtain 
{wn+i + r„ - d„_i)+ w„. □ 

Proposition 4.6. Let the sequences of r.v.'s (a„)„gN-> and (s„)„gN* be as in j j^. J| or as in The 

conditional law of (X]"=i ^i; X]r=i ■^i)' given that {X]i=i ^ — ^^1 the same as the 

unconditional law of (maxi<j<„{^^^^ + X]"=^/+i iiiini<i<n{I]i=2 '^^ + J27=j+i '^^})■ 

Proof. Set ti„ ~ Wn + Sji = Wn+i + (the sojourn time of customer n in the queue). By developing 
(fTS)) . we obtain 

fc fc-i 

_i +r„ = sup{^ri -^di} . (17) 



By summing the equalities (|16p . we obtain on the event {wi = 0}: 



i=l 1=1 
n 

= Vaj+ max{ V - V aj = max{Vaj+ V sj 

z — l<')<n -"^^ — ' -"^^ — ' l<7<n ^ — ^ ^ — ^ 



n n+1 n j n+1 



I<j<n ^ — ^ ^ — ^ l<jl 

and summing Q on the event {t;i — 0}, we get: 

n n 
y^^n = ^Si- Wn+1 
i=l i=2 

n n n j n 

ESj - max { Si - fli} = min {Y^ +^"04}. 

l<i<n l<i<ri '^-^ 

4=2 ~ i=i+l j=i+l ~ i=2 i=j+l 

To summarize, we have 

n j n+1 n J n 

" S + H XI " ,3^^ S ■^^ + X! "^j' ■ (^^) 

— ' l<j<n ^ — ' ^ — ^ ^ — ^ ISjS^ — — 

i—1 i—1 i—j + ^ i—1 i—2 i—j+1 

Applying Theorem l4.1l or H^ the distribution of (X]"=i '^j; J27=i given that {J2'i=i — J2i=i '^^ — 
0} is the same as the law of '^O given that {J2'i=i d-i — J2i=i '^ij ^ 0}- By (fTTl) . the 

event {J2i=i — J2'i=i Vfc > 0} is equal to {vi — 0}. 

Using p8|) . we have that the distribution of {J27=i^i^ S"=i ^i) given that {vi = 0} is the 

same as the law of (maxi<j<„{X]Li + I]r=i+i ■^O; n^iiii<i<"{Si=2 + Z^iLj+i " J) given that 

{vi — 0}. This last conditional law is the same as the unconditional law of (maxi<j<„{^^^j^ + 

J2'i=j+i Si},mini<j<n{I]i=2 +I]"=j+i "i})- Indeed, the variables {ai)i>i and {si)i>2 are independent 
of = 0}. 

By fitting all the parts together, we obtain the desired result. □ 

It is possible to extend Proposition l4.6l to the limiting case E[ai] — E[si], and also to higher dimensions, 
by adapting the methods of [25l [181 HI] to the present setting. 
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5 Transient Behavior and RSK Representation 



In Sections 15.11 and 15. 3[ the model described and the results obtained have a purely combinatorial 
nature. In particular, the results are true pathwise, without any probabilistic assumption. These results 
are then used in Section 15.41 for geometric queues or stores: here the pathwise arguments are coupled 
with probabilistic ones. 



5.1 The saturated tandem 

We consider another aspect of the dynamic of queues and stores: the transient evolution for the model 
starting empty. More precisely, consider the model of fJH] under the assumption that wq = Aq ^ sq ^ 
(which implies Dq = tq = 0) and focus on the customers, resp. time slots, from 1 onwards. 
It is convenient to describe such a model with a different perspective. We first do it for the queue. 
View the arrivals as being the departures from a virtual queue having at instant an infinite number of 
customers (labelled by N*) in its buffer. The service time of customer n in the virtual queue is a„_i. 
We describe this as a saturated tandem of two queues. 

Let us turn our attention to the store. View the supplies as being the departures from a virtual store 
having an infinite stock at the end of time slot 0. In the virtual store, the request (=departure) at time 
slot n is s„ . This is a saturated tandem of two stores. 

Now we want to fit these two descriptions together. Denote the virtual queue/store as queue/store 1 
and the other one as queue/store 2. For convenience, set u{n, 1) = a„_i and u(n, 2) = s„ for all n > 1. 
The saturated tandem is completely specified by the family (u(n, j), n g N*, z = 1, 2) of input variables. 
These variables are the services, resp. requests, when the model is seen as a tandem of queues, resp. 
stores. Next table gives the input variables in the saturated tandems of two queues/stores. 



Customer / Time slot 


1 


2 




n 


Queue 2 / (Virtual) Store 1 


u(l,2) = si 


u(2,2) = S2 




u{n, 2) = Sn 


(Virtual) Queue 1 / Store 2 


u(l, 1) = ao 


u(2,l) = ai 




u{n, 1) = a„-i 



A couple of observations are in order. Observe that u(n, i) is the service of customer n in queue i, and 
the request at time slot n in store 3 — i. In other words, the elements (queues or stores) associated with 
a given sequence (m(-, i)), i = 1,2, are crossed in reverse orders in the queueing/storage tandem. This is 
illustrated in Figure [4] (set K = 2). Observe also that there is a shift in the time slots for the storage 
model: the departure from store 1 at time slot n is the supply of store 2 at time slot n+1 (contrast this 
with the situation for the queues). This is consistent with a model in which we view a time slot as being 
decomposable in three consecutive stages: first, the supply arrives; second, the request is made; third, 
the departure occurs. See Figure [2l 



Queue 1 ^ Queue 2 

u(.,l) u(.,2) 



Queue K 
u(.,K) 



Store K 



Store K-1 



Store 1 



Figure 4: Queues and stores are gone through in opposite orders 

The above setting is naturally extended to define the saturated tandem of K queues/stores. Such a model 

is entirely defined by a family of N- valued variables (m(^, j))ieN*,ie{i k}- (In 33 the variables were 

valued in N*. This restriction is not necessary here.) For the queueing model: (i) at instant 0, queue 
1 has an infinite number of customers labelled by N* in its buffer, and the other queues are empty; (ii) 
u{i,j) is the service of customer i at queue j; (iii) the instant of departure of customer n from queue i is 
the instant of arrival of customer n in queue i + 1. For the storage model: (i) at the end of time slot 0, 
store 1 has an infinite stock and the other stores have an empty stock; (ii) u{i, K + I — j) is the request 
at time slot i in store j; (iii) the departure at time slot n from store i is the supply at time slot n + 1 in 
store i + 1. The models are depicted in Figure [H 



10 



5.2 Interacting particles representation 



It is helpful to describe the above tandem models as interacting particle models (see for instance T9]). 
Commonly, interacting particle systems are defined as continuous-time Markov processes. In contrast, 
the description below is pathwise and the time is discrete. 

Classically, the queues in tandem may be viewed as the following zero-range process. The set of sites is 
{!,..., -ftT} and a site is occupied by a number of particles in N U {+00}. At each site, the front particle 
(if any) has a clock whose value is decremented by 1 at each time slot. When it reaches 0, the particle 
jumps to the site on its right. The right-most site has an infinite number of particles. This is illustrated 
in Figure [5l The variables {u{i,j))i correspond to the initial values of the clocks of the successive front 
particles at site j. 



Instant (T+1)- 



Figure 5: The zero-range process for queues in tandem 

We now describe the stores in tandem as an original zero-range model: the bus-stop model. At the 
beginning of each time slot, a bus arrives at each site. The size of the bus varies from site to site, and 
from time slot to time slot. The particles take place in the bus until it is filled in, and are transported 
to the next site on the left, during the current time slot. The left-most site has an infinite number 
of particles. This mechanism is illustrated in Figure [6l Here the variables {u{i,j))i correspond to the 
successive bus sizes at site j. 



Figure 6: The zero-range process for stores in tandem 



Alternatively, the above two models can also be viewed as exclusion processes. The set of sites is 
embedded in N and each site now contains either or 1 particle. 
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The queues in tandem are associated with the classical TASEP or totally asymmetric exclusion process. 
For the stores in tandem, the exclusion process is as follows. Particles take jumps of integer length to the 
right which are constrained by the interdiction to overpass other particles. In one time step, particles 
jump in order from left to right. This model is close but different from the exclusion processes considered 
in IMIBSI. 



Queues 



Stores 



\ \ 



• • o o • 



• o o o • • 



Figure 7: Exclusion processes for the queues and the stores 



The mapping between zero-range process and exclusion process, both for the queues and the stores, is 
illustrated in Figure [71 



5.3 Robinson-Schensted-Knuth representation 

A partition of n G N* is a sequence of integers A — (Ai, . . . , Afe) such that Ai > • • ■ > Afc > and 
Ai + • • ■ + Afe =71. We use the notation A h n. By convention, we identify partitions having the same 
non-zero components. The (Ferrers) diagram of (Ai, . . . , Afc) h n is a collection of n boxes arranged in 
left-justified rows, the i-th row starting from the top consisting of Ai boxes. A (semi-standard Young) 
tableau on the alphabet {1, ...,£} is a diagram in which each box is filled in by a label from {1, ...,£} in 
such a way that the entries are weakly increasing from left to right along the rows and strictly increasing 
down the columns. The shape of a diagram or tableau is the underlying partition. A standard tableau 
of size n is a tableau of shape A h rt whose entries are from {1, . . . ,n} and are distinct. In Figure [8] we 
have represented a diagram on the left and a tableau of the same shape on the right. 



112 3 

2 2 

3 4 
4 



Figure 8: Ferrers diagram and semi-standard Young tableau of shape (4, 2, 2, 1) h 9 

The Robinson-Schensted-Knuth row-insertion algorithm (RSK algorithm) takes a tableau T and i E N* 
and constructs a new tableau T i. The tableau T ^ i has one more box than T and is constructed as 
follows. If i is at least as large as the labels of the first (upper) row of T, add a box labelled i to the end 
of the first row of T and stop the procedure. Otherwise, find the leftmost entry in the first row which is 
strictly larger than i, relabel the corresponding box by i and apply the same procedure recursively to the 
second row and to the bumped label. By convention, an empty row has label 0. With this convention, 
the above procedure stops. 

Consider a word v — vi ■ ■ ■ Vn over the alphabet {1, . . . , k}. The tableau associated with v is by definition 

P = (• • • ((To ^ Vi) ^ «2) ■ • • ^ Vn) , 

where To is the empty tableau. Observe that P has at most k non-empty rows. Classically, the length 
of the top (and longest) row of P is equal to the longest weakly-increasing subsequence in v. 

Remark 5.1. While building the tableau P, it is possible to build another tableau of the same shape 
in which the entries, labelled from 1 to |t;|, record the order in which the boxes are added. This 
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recording tableau is a standard tableau of size \v\. By doing this, one defines a bijection between words 
of {1, . . . , fc}" and ordered pairs of tableaux of the same shape, the first being semi-standard over the 
alphabet {1, . . . , fc} and the second being standard and of size n. This bijection is often referred to as 
the Robinson-Schensted-Knuth correspondence. Here, we do not need this result and we do not consider 
the recording tableau. 

Consider a family U = {u{i, j)){ij)^\i^,,,^M}K{i....,K} of variables valued in N. (Here we do not make any 
assumption on the variables, which may or may not be random.) We associate with U the word w(U) 
over the alphabet {1, . . . , K} defined by 



w{U) = wi . . . Wn, with Wi = ^J^^^ ' ■ ' ' J 

u{i, 1) u{i, 2) • • • u(i, K) 



(19) 



Set M = \w{U)\ = X;^=l Ejli w(i, j)- For i = 1, . . . , X, define 

Vn < M, Xi{n) — \{i<n \ w{U)j = i\\ . 

Given two maps x,y : {1, . . . , N} — > N, define the maps xVy, x A y : {1, . . . , N} — > N as follows: 

Vn < N, 



x\/y{n) — max x(m) + y{n) — y{m) , x A y{n) — min x{m) + y{n) — y{m) 

0<m<n L J 0<m<n 



(20) 



Denote by P{U) the tableau obtained from w{U) by applying the RSK algorithm and let (Ai, . . . , A^) 
be its shape. The following holds 

Ai = a;iVx2V . . . Va;if (Af), A^ = A . . . A 2:2 A a;i(M) , (21) 

where the operations V, A are performed from left to right (the operations V, A are non-associative). 
The expression for Ai follows from the fact that Ai is the longest weakly-increasing subsequence in w(U). 
The expression for Xx is proved in 2^ Theorem 3.1]. In fact, the result from |23j is more general: there 
exists a min- max- type operator Tk such that (Ai, . . . , Aa') = Tk{xi, . . . , xk)- 

A lattice path is a sequence tt ~ ((*i,ii), ■ • ■ , with {ik,jk) G The steps of tt are the differences 

(ik+i - ik,jk+i - jk), fc = 1, . . . / - 1. 




Let H be the set of lattice paths from (1,1) to {N, K) with steps of the type (1, 0) or (0, 1). All the paths 
in H have the same number of nodes: N + K —1. Let H be the set of lattice paths from one of the points 
in {(1 -f- i, if — i), i = 0, . . . , if — 1} to one of the points in {{N — j,l + j), j — 0, . . . , N — 1} with steps 
of the type (1 -I- i, —i), i >0. 
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Observe that all the paths in 11 also have the same number of nodes: {N ~ K + 1)+. In particular, 11 is 
empty when N < K. An example of a path from each of the two sets is provided in Figure [9l 

The following result holds. When interpreting it, recall that queues and stores are arranged in opposite 
orders in the tandems, see Figure [H It is also fruitful to compare the statements of Proposition 14.61 and 
Theorem O for K = 2. 

Theorem 5.2. Consider a saturated tandem of K queues/stores with variables u{i,j),i G N*,j G 
{!,... ,K}, as defined in ^5.1[ Fix iV e N*. Let D be the instant of departure of customer N from 
queue K and let R be the cumulative departures over the time slots 1 to N in store K. Set U = 
{u{i, j))i<ci<N , i<j<K o,nd define w{U) as in \19\) . Let (Ai, . . . , A^) be the shape of the tableau associated 
with w(U). We have 



Ai = max u{i,j] 



min V u{i,i 



D 



R. 



(22) 
(23) 



First, we state and prove a preliminary result. 
Lemma 5.3. We have 



xk a ... a X2 a Xi{M) = min > u{i, 



(24) 



Proof. By definition. 



Xk A . . . A xi{M) — min xxirriK) + XK-i{mK-i) - XK-i{mK) + ■ 

l<mK'<---<mi — M 



xi{mi) ~ xi{m2) 



Let (mjf , . . . , mj) be a sequence of integers realizing the minimum in the above equation. Consider the 
set of integers 

k A'-l 

I^{y^w^\ + y. <k + 1, j), e {0, . . . , X - 1} . 



i=l 



We may assume wlog that m*j^ E I. Indeed assume that mj^ ^ /, and let uk be the smallest integer 
such that hk G L^uk > m}^- Then, we clearly have xxinK) — XK{'m*j^) and 

Xk A . . . A xi{M) - XKin-K) < min XK-iiniK-i) ~ XK-iinx) ^ h a;i(M) - Xi(m2) 

riK <mK-i<---<m2 

< min XK-i{mK-i) ~ XK-i{m*K) + ■ ■ ■ + Xi{M) ~ Xi{m2) 

= Xk A . . . A xi{M) - XKi'm}^) . 



So we assume that £ /. Let iK be such that to|^ = ^ l^d + X^jLi^ j)- We have 

XK(m^) = ^2]=!^ "^{.h ^) ■ Then we prove with the same type of argument that we can assume wlog 
thatm|f_i e {ELi\^^\ + Ef=i^uik + l,j), ke{iK,---,K-l}}. 



w(U) : 



. (K-1)(K-I)(K-1)KKK ...Kill KKK 



Figure 10: Illustration of the proof of Lemma 15.31 

Since m'^ G /, we have XK{Tn*K) — "(^j j) for some iK- Then we prove with the same type of 

argument that we can assume wlog that m*^_^ G {X]i=i l^d + SjL^^ u{k,i), k G {iK + 1, . . . , if — 1}. 

It implies that XK-iirriK-i) — XK-iiiriK) — Sj^^A+i ""(-^ ^ ^^j) fo'" some iK-i- By repeating the 
argument, we obtain the desired result. □ 
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Proof of Theorem \5.2l The result for Ai is essentially due to Schensted. The left equality in ((22|) is ([21]) 
and the right equality appears for instance in [211 123 E] • As far as Xk is concerned, Lemma 15.31 leads 
to 

= a; A' A ... A 2:2 A xi{M) = min >^ . 

TrenL."^ -I 

It remains to prove the right equality in (|23p. Let us denote by r(i — the amount departing at time 
slot i from store j and w{i,j) the quantity remaining at the end of time slot i in store j. According to 
©, we have 

r{n, k) ~ r{n — 1, fc — 1) + w{n, k) — w{n + 1, fc) , (25) 
Applying ^ to our tandem model, 

w(n+l,fc) = [w{n,k) + r{n - l,k - \) - u{n,K + 1 - k)\^ 

n — 1 71 

= max [^(r(i,A; - 1) - ^ X + 1 - fc))] . 

~ ~ i—m 2— m+1 

Using ([25]) . for n > iiT, we get 



r{n,K) = r(n - 1,X - 1) + u;(n,ii:) - + 1,X) 

= r(n - 2, X - 2) + w(n - 1, if - 1) - w(n, K - 1) + w{n, K) - w{n + 1, K) 

= r(7i - if + 1, 1) + ^ [w{n -i,K -i) - w{n - i + \,K - i)] 

i=0 

= u{n - if + 1, if ) + ^ [w{n ~ i,K - i) - w{n - i + \,K - i)] . 

i=0 

Hence, 

N N 

i?=^r(n,if) = 51 ^("'^) 

n=l n— 

AT N K-1 

= ^ M(n - if + 1, if ) + ^ ^ w(n -i,K -i) - w{n - i + 1, if - z) 

N-K+1 K-1 N 

= if) + 5^ w{n — i, K — i) — w(n — i + 1, if — z) 

n— 1 i—0 n—K 

N-K+1 K-1 

u{n,K)+Yw{K-i,K-i)-w{N-i + l,K^i). 

n=l i=0 

Now recall that w{k,k) = 0, for all fc > 1. We obtain 

N-K+1 K K 

R = ^ u{n,K) + Yw{i,i) ~Yw{N - K + i + 

n—l i—1 i—1 

N-K+1 K 

^ u(n,K)^Yw{N ~ K + i + . (26) 

n— 1 i—1 

Moreover, set A — w{N + 1, if) + w{N, if — 1), we have 
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N-1 N 



A = ^ma^^^^ r{n, K - I) - ^ u{n,l) +w{N,K-l) 

^ n—ii n— + 1 

N-1 N 

= ^max^^^[r{n -1,K ~2) +w{n,K - I) -w{n + l,K -1)]- ^ u{n,l) +w{N,K-l) 

^ n—ii n— + 1 

N-2 N 

= ^max^[ r{n,K -2) + w{ii,K -I) -w{N,K -I) - ^ u{n,l) +w{N,K~l) 

^ n—ii—l n—ii+1 

N-2 N 

= max^l^ r{n, K — 2) + w{ii, K — I) — u{n,l) 

^ n—ii — l n— 

N-2 ii-2 N 

= max I r(n,K-2)+ max r(n, K - 2) - u(n,2) - u(nA )] 

n—ii — l n—L2 71—12 + 1 n— 

11—12 11—12 + 1 n— 

Applying this repeatedly leads to 

K N-K+l iK-2-l N 

Yw{N - K + i + = max^^l^ ^ u{n,K)- ^ u{n,K-l)^ h ^ M(n, 1) . 

-i— 1 — -ff 1 1— n—i-K — i n— iit-i+l n— ii + l 

(27) 

Combining (pS)) and (|?7|) . we obtain 

jjf-l — 1 i/<-2 — 1 iV 

i? = min I" "V ui'n,K)+ \^ u(n, K - 1) ^ h u(n,l) 

n— 1 n— i/c — 1 + 1 n— H+l 

□ 

Let us comment on Theorem 15. 21 First, it is interesting to observe that the 6 quantities appearing in the 
equalities ([^ and (^5)) all come from the same model. Now, consider the two identities ((^ and ((^ 
separately. The first one is classical and has led to several interesting developments for queues in series. 
A by-product of the newly proved identity (|23|) is therefore to pave the way to obtain analog results for 
stores in tandem. We set out this program in some details in the next Section. 

5.4 Stochastic models, interchangeability, and hydro dynamic limits 

If r is a tableau over the alphabet {!,..., /c}, we write — JliLi^rN where ai is the number of 
occurrences of the label i in the tableau. The Schur Junction s\ associated with the partition A = 
(Ai, . . . , Afc) is defined by 

Sx{Xi,...,Xk) = ^ x"^ , 
T: sh(T) = A 

where sh(T) is the shape of the tableau T. We refer the reader to the books [321 131] for rnore about 
Schur functions and their connection to the RSK algorithm. 

Suppose that the random variables U = j))(ij)e{i,...,jv}x{i,..../<'} ^-re independent and that u{i,j) 
is geometrically distributed with parameter q^, for some fixed {qi,...,qK) e (0, 1)^. Then ^ the 
law of the random partition A = (Ai, . . . , A^) associated with w{U) is given by 

P{A = /} = a(g)^sz(g)si(l,...,l), I Pk , (28) 
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where Pk is the set of mteger partitions with at most K non-zero parts, a{q) — 11^(1 ~ 1j) ^-i^d si is the 
Schur function associated with the integer partition I. In particular, the law of A is symmetric in the 
parameters qi, . . . ,qK- 

This extends to the level of processes: if we write A — A^^-' to emphasize its dependence on N, it is easy 
to check [52] that the random sequence is a Markov chain in Px with transition probabilities given 

by 

P{A(^+i)=Z| AW=m} = a(g)^, (29) 

Srn{q) 

for all m and I such that h > mi > /2 > rn2 > ■ ■ ■ . In particular, we can see that the law of the sequence 

is symmetric in the parameters qi,. . . ,qK- 
These remarks give the law of the departure processes from both of the stochastic tandem queueing and 
storage models described above. Consider a saturated tandem of K queues/stores as defined in i j5.ll 
Assume that the random variables u{i,j), i G N*,j G {1, . . . , K} are all independent, and that u{i,j) is 
geometric of parameter qj G (0, 1). Let D'-^^ be the instant of departure of customer N from queue K 
and let i?^^) be the cumulative departures over the time slots 1 to in store K. By Theorem l5.21 D'-^^ 
and i?*^^) are, respectively, the lengths of the longest and shortest rows A^^"* and A^'' of the random 
tableaux obtained by the RSK algorithm. 

Interchangeability. A result by Weber [38] (see also Anantharam [1]) states that a series of -/M/l 
queues arranged in tandem are interchangeable, i.e. the output process has the same distribution under 
any reordering of the queues. This is exactly the fact that the law of D is symmetric in the parameters 
and can be seen now as a direct consequence of the symmetry of the Schur polynomials. We can now 
extend this result as follows: 

Theorem 5.4. The law of {D'-^\ R'-^^)n is unchanged if the parameters of the geometric laws are 
modified from (gi, . . . , qx) to (ga-(i)i • • ■ j Q.a(K))j for any permutation a of {1, . . . , AT}. 

In particular, stores arranged in tandem are interchangeable. 

Hydrodynamic limits. Proving 'shape theorems' (also known as 'hydrodynamic limits') and fluctu- 
ation results is an important issue in interacting particle systems, see for instance [16]. The law of D^^^ 
is well-known and the connection with RSK has been exploited before to obtain detailed asymptotic and 
fluctuation results for the corresponding interacting particle system. The tandem storage (and bus-stop) 
model however has not been discussed in the literature. A by product of our analysis is that this model 
is amenable to some very precise analysis and development. 

Indeed, consider Theorem l5.4l and the important special case when the parameters are equal, qj = q say. 
Then the random variables D and R arc jointly distributed as the largest and smallest 'eigenvalues' of 
the Meixner ensemble. In the exponential case, replace Meixner by Laguerre. For example, if the u{i,j) 
are exponentially distributed with parameter 1, then the cumulative departures i?^^^ over the time slots 
1 to N in the corresponding tandem of K stores is distributed as the smallest eigenvalue of the K x K 
random matrix AA* , where A is a, K x N matrix with i.i.d. standard complex normal entries. This 
distribution is known, see for example [29| . When N — K it is exponential with mean K. This can be 
used as a starting point for the analysis of the fluctuations and hydrodynamic limit corresponding to the 
bus-stop model of Section [5.21 a programme which is beyond the scope of this paper. 

Fixed point results. Finally we remark that the identity [ min^^jj ^-^g^ u(i, j) = i?] can also be 
used as a cornerstone for proving fixed point and related weak convergence results for stores in tandem. 
In the context of queues in tandem, these questions have been pursued for some time and the identity 
[maxjren j)e7r = ^] ^ ingredient in solving them [51 \7[ \^\T7\ . 
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